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Abstract. The Kramers-Kronig relations are derived for the permittivity of the 
usual plasma model which neglects dissipation and of a generalized model which takes 
into account the interband transitions. The generalized plasma model is shown to be 
C*~) ' consistent with all precision experiments on the measurement of the Casimir force. 
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^ ' PACS numbers: 05.30.-d, 77.22. Ch, 12.20.Ds 

Qh, 1. Introduction 

+-> 

It is common knowledge that the Kramers-Kronig relations connect real and imaginary 
parts of an analytic function describing some causal physical process. In statistical 
physics and electrodynamics any material susceptibility satisfies the Kramers-Kronig 



relations [H Specifically, they are used to calculate the real part of the dielectric 
permittivity e'(u), along the real frequency axis and the dielectric permittivity, e(i£), 
along the imaginary frequency axis [2j [3]. Both are expressed through the imaginary 
part of the permittivity, s"(u), at all real frequencies. 

In the last few years the Kramers-Kronig relations have been repeatedly used to 
calculate the thermal Casimir force in the framework of the Lifshitz theory (see, e.g., 
reviews [HE] and recent proceedings [6]). The Casimir force [7] acts between neutral 
material bodies and originates from the zero-point oscillations of the electromagnetic 
field. The Lifshitz theory allows one to express the Casimir force at a temperature T in 
terms of e(i£) of the body materials at discrete Matsubara frequencies £j = 2nk B Tl/h, 
where ks is the Boltzmann constant and / = 0, 1, 2, . . . . The tabulated optical data 
for the complex index of refraction and hence for e"(uj) are, however, available only in 
some restricted frequency regions [8]. Therefore the direct calculation of e(i£) using 
the Kramers-Kronig relations is not possible and different approaches to find e{iC,) have 
been proposed. 

In the first approach [9l [10] the quantity e"(uj) obtained from the tabulated optical 
data is extrapolated using the imaginary part of the Drude dielectric function to all lower 
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frequencies including zero frequency. Then e(z£j) at all / > is found from the Kramers- 
Kronig relations and is substituted into the Lifshitz formula at nonzero temperature. 
From the theoretical point of view this approach may seem straightforward, but it leads 
to a violation of the Nernst heat theorem for perfect crystal lattices with no impurities 
[IT] and to contradiction with the experiment measuring the Casimir force at separations 
from 160 to 750 nm |12j . The second approach [T2l [T3] is based on the concept of 
the Leontovich surface impedance [2]. This approach leads to practically the same 
contributions to the Casimir force, as the first approach, at all Matsubara frequencies 
with I > 1. The contribution of the zero-frequency term is, however, different and fixed 
by the impedance used. As a consequence the second approach is not applicable at 
separations below the plasma wavelength (equal to 137 nm for Au) where the Leontovich 
impedance boundary condition becomes invalid. The third approach [HI [15] does not use 
the tabulated optical data but employs the dielectric permittivity of the free electron 
plasma model at all frequencies. Both the second and the third approaches are in 
agreement with thermodynamics. They are also consistent with the experiment [12] 
performed at separations above the plasma wavelength. However, both the second and 
third approaches cannot be applied in the experiment [JU [17] where measurements start 
at short separations of 60 nm. This experiment although performed at T = 300 K was 
found to be consistent with the Lifshitz theory at zero temperature. (The comparative 
analysis of all approaches is contained in [18-20].) Note that the third, plasma model 
approach, may seem to be in disagreement with the Kramers-Kronig relations because 
the dielectric permittivity of the plasma model is entirely real. In this connection the 
plasma model approach has been criticized [20] for the complete neglect of dissipation. 
Thus at the moment none of the theoretical approaches to the thermal Casimir force is 
consistent with all the available experimental information. 

In the present paper we derive the generalized Kramers-Kronig relations for the 
permittivities of the free electron plasma and a plasma-like model which incorporates 
dissipation due to interband transitions. We demonstrate that the permittivity of the 
plasma model (as any function analytic in the upper half-plane) satisfies the Kramers- 
Kronig relations if the contribution from the pole of the second order at zero frequency is 
correctly taken into account. Then we compare theoretical computations of the thermal 
Casimir force using the free-electron plasma model and the generalized plasma model 
incorporating interband transitions with the zero-temperature Casimir force calculated 
using the tabulated optical data. We demonstrate that the theoretical results using the 
generalized plasma model are in good agreement with experiment. Thus currently it is 
the only model for the thermal Casimir force which is consistent with all measurements 
performed to date. We conclude with a discussion of different types of dissipation 
processes and their role in the theoretical description of the Casimir force. 
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2. Kramers-Kronig relations for plasma and plasma-like models 

We consider the generalized plasma-like dielectric permittivity of the form [21 



e(u) - 1 = A(u) - j£ (1) 



where u p is the plasma frequency and the oscillator term 

A M = £ U 2 j wu ( 2 ) 
j=l U j ~ u ~ l 9j UJ 

takes into account the interband transitions of core electrons. Here Uj ^ are the 
resonant frequencies of the core electrons, gj are the respective relaxation frequencies, 
fj are the oscillator strengths, and K is the number of oscillators. The dielectric 
permittivity (0Q), (j2J) was used in Sec. 7.5(D) of [21] for the description of a metal 
at frequencies much larger than the Drude relaxation frequency. The term —uj 2 /uo 2 in 
Eq. ([1]) describes the free conduction electrons and leads to a purely imaginary current. 
This contribution to e{u) is entirely real and does not include dissipation. It must 
be emphasized that the oscillator term ([2]) does not include the oscillator with zero 
resonant frequency Uq — 0. Thus it does not describe conduction electrons but only the 
core electrons. If the core electrons were excluded from our consideration then fj = 0, 
A(uj) = 0, and the dielectric permittivity (CD) leads to the usual plasma model. Note that 
for the purpose of the computations below we follow the notations from [22] (Level 2, D) 
for the parameters of the interband oscillators. Because of this we have replaced the 
relaxation parameter Tj in [21] for gj and the oscillator strengths AirNe 2 fj /m, where N 
is the number of molecules per unit volume, as in [21], for fj. Here we also use 1 + A(u) 
in place of £&(o>) [21]. Equations ([I]) and (J2J) incorporate dissipation due to interband 
transitions but do not include processes of electron scattering on phonons, impurities, 
grain boundaries, surfaces and other electrons. Below we investigate the mathematical 
properties of Eqs. (0Q), (j2J) for the complete frequency range from zero to infinity. The 
physical justification for the choice of e{uS) in Eq. p]) is discussed in Sec. 4. 

The characteristic feature of the dielectric permittivity (OQ) is the second order pole 
at zero frequency. Let us demonstrate that ([T|) satisfies the Kramers-Kronig relations in 




Figure 1. The integration contours (a) in Eq. ([3]) and (b) in Eq. (fTTj) consisting of 
the real frequency axis and the semicircle of infinitely large radius. 
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both cases A(uj) = and A(uj) 7^ 0. For this purpose we consider the integral 

/ e M^±du, = 0i (3) 

Jd UJ — UJq 

where ujq is real and the contour C\ is presented in Fig. 1(a). Inside G\ the function 
under the integral is analytic and thus the equality d3J) follows from the Cauchy theorem. 
At infinity e{ui) — > 1 and the function [s(uj) — l]/(u> — uj ) therefore tends to zero more 
rapidly than 1/uj. Because of this the integral along the semicircle of infinite radius is 
zero. We pass around the points and ujq along the semicircles C p and C$ with radii p 
and S, respectively. It is easily seen that 



e{uj) 



-duj 



-7uRes- 



eM-1 



IC S UJ — UJq UJ — UJq 

The similar integral around the point is more involved. Using Eq. ([1 
as a sum of the integral 

A(uj) 2,4(0) 



m [e(uj ) - 1] . (4) 
we represent it 



-duo 

IC P LO - UJq 

which vanishes when p — > 0, and 



-p, 



C p U 2 (u - Uq) Uq Jc p 



UJQ 

UJ 2 



OJ 



UJQ 



1 

UJ 



du. 



(5) 



(6) 



Direct integration along the semicircle C p results in 



du 



2 

: P> 

C p UJ - LUq UJq 

duj 2 

— = — = UJ P 
uj z p 



du 

( ■„ UJ 

duj 



-7TZ, 



-00 UJ Z [UJ — UJq) 

where the integral is taken as a principal value. [Note that the last integral cannot be 



evaluated as in Eq. (j4j) because Res(l/c<j z 



u=0 



and both integrals around the upper 



and lower semicircles are divergent and opposite in sign.] 
Substituting Eqs. (|I])-([7j) in Eq. Q we arrive at 



ITXUJl 



V 



UJt\ 



i] +P J_ 



duj 



—00 UJ — UJq 



e\uj 



-1 + ^ 

UJ 1 



0. 



(8) 



Now we replace the integration variable uj by £, ujq by uj, and represent the function 
e{yj) in the form of e(u) = e'(uj) + ie"{uj). Taking into account that 

duj 



P 







-00 UJ — UJq 

and separating the real and imaginary parts in Eq. 
Kramers-Kronig relations 



s\uj) = 1 + -P 



7T 



Z-uj 



d£ 



e"(u) 



(9) 

we obtain the generalized 
e'(0 + 



-P 



7T 



e 2 



Note that the standard relations [2] obtained for permittivities with no pole at uj 
do not contain terms uJp/C, 2 on the right-hand sides of Eqs. (fTDl . 



d£.(10) 




Kramers- Kronig relations for plasma-like permittivities 



5 



The dielectric permittivity along the imaginary frequency axis can be determined 
through the use of the integral 
r u[e(u)-l] 
Jc 2 



-duj 



m [euu>Q) 



1] 



UJ 2 + uj 2 

along the contour C2 in Fig. 1(b). By integrating over C2 we get 

uj [e(u) - 1] 



UJq J-. 



-did = in [e(iuJo) — 1] . 



(12) 



UJ 1 + UJq 

Now we make the same replacement of variables as above, separate the real and 
imaginary parts of e{uS) under the integral and use the identities 



P 



0. 



p 



50 fr'(0 

oo£ 2 +U 2 



di = 0. 



The result is 



e(iuj) — 1 



-P 



TT 



00 ^ 2 + U) 



UJ 



(13) 



(14) 



For the usual plasma model e"(u) = 0, e'(uS) = 1 — u> 2 /u> 2 , e{iuS) = 1 + u 2 /u 2 , the 
generalized Kramers-Kronig relations (fTUj) . (fT4l) are satisfied with the use of Eq. (Q. On 
the contrary the same plasma model violates the standard Kramers-Kronig relations. 
Note that sometimes [23] the plasma model is ascribed a nonzero imaginary part 



e"(u) 



-— lim 



%5H 

UJ 



(15) 

uj g~*o uj + ig uj 

which is obtained from the Drude dielectric function in the limit of zero relaxation 
parameter. This makes it possible to formally satisfy the standard dispersion relation 
for the dielectric permittivity along the imaginary frequency axis j2] given by Eq. f|T4l) 
without the uj 2 /uj 2 term. However the other two standard Kramers-Kronig relations 
with e"(uj), as given by Eq. ( 1T51) . become meaningless. Thus the permittivity of the 
collisionless free electron gas is entirely real. 

It is easily seen that the plasma-like dielectric permittivity (CD), (J2J) satisfies the 
generalized Kramers-Kronig relations ffTUj) and f|T4l . This can be verified by direct 
substitution. For example, the substitution of Eqs. (0Q) and (j2D in the first equation of 
ffTUj) leads to 



K 

E 



UJ' 



1 



K 



j=l [UJ 



UJ' 



K 



g]uj 2 
fj9j 



UJ) 



7T 



(-1 



UJ' 



UJ > 



+ gy 

dy 



(16) 



»2/ 4 - 2/3 jV 2 + 1 



UJi 



y 2 dy 



where (3j = 1 — g 2 /(2uj 2 ). When the following: 



»2/ 4 - 2P 2 + 1 



7T 



00 2/ 4 - 2^-y 2 + 1 J-00 - 2/3 j2 / 2 + 1 J2(l-^ 



;i7) 
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is taken into account Eq. (fTBl) is satisfied. The second equation in (fTU|) and Eq. (fT4"|) can 
be verified in a similar way. 

3. Calculation of the Casimir force using the generalized plasma model and 
comparison with experiment 

We have calculated the thermal Casimir force F th acting between Au coated test bodies 
in the most precise short-separation experiment [16J (a sphere of R = 95.65 fim radius 
and a plate) by substituting the generalized plasma dielectric function ([[]), ([2]) into the 
Lifshitz formula. The Lifshitz formula for the force between a sphere and a plate was 
obtained by means of the proximity force approximation (PFA) as 2nRJ-' where T is 
the free energy in the configuration of two parallel plates (see [5], [17] for details). Note 
that recently the high accuracy of PFA was confirmed using the path-integral approach 
[24"1 [25] and worldline numerics [26] . The oscillator parameters for Au in Eq. ([2]) were 
found in [22j [27] using DESY data. They are as follows: K = 3, uj% — 3.87 eV, 
fx = 59.61 (eV) 2 , gi = 2.62 eV, to 2 = 8.37 eV, f 2 = 122.55 (eV) 2 , g 2 = 6.41 eV, 
co>3 = 23.46 eV, f 3 = 1031.19 (eV) 2 , g 3 = 27.57eV. Computations were performed at 
the laboratory temperature T = 300 K at different experimental separations a with 
Up = 9.0 eV [8j [T7J. The obtained magnitudes of the Casimir force are presented in 
Table 1 (column 2). For comparison in column 3 we present the force magnitudes 
obtained using the usual plasma model of Eq. (0Q) with A{u) = 0. Column 4 lists the 
force magnitudes from the zero-temperature Lifshitz formula using the tabulated optical 
data for the complex index of refraction (recall that in [161 E] the experimental data 
were compared with theory at zero temperature). As is seen from Table 1, at short 

Table 1. Magnitudes of the Casimir force at different separations in column 1 
computed at T = 300 K using the generalized plasma- like model (column 2) , the usual 
plasma model (column 3), and at zero temperature using the tabulated optical data 
for the complex index of refraction (column 4). 



Force magnitude (pN) 



ft 


Generalized 


Plasma 


Force at zero 


(nm) 


plasma model 


model 


temperature 


60 


531.1 


483.2 


527.4 


70 


358.8 


332.2 


356.1 


80 


254.9 


239.1 


252.8 


90 


188.2 


178.3 


186.5 


100 


143.3 


136.8 


141.9 


120 


88.94 


86.00 


88.01 


150 


49.30 


48.19 


48.71 


200 


22.75 


22.46 


22.44 


250 


12.37 


12.28 


12.19 


300 


7.478 


7.438 


7.355 
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Pf - (F^>) (pN) F; h - (F exp ) (pN) 




100 150 200 250 300 100 150 200 250 300 



z (nm) z (nm) 

(a) (b) 
Figure 2. Differences between the theoretical and mean experimental [16] Casimir 
forces versus separation. Theoretical forces are computed using (a) the generalized 
plasma-like model and (b) the usual plasma model. 



separations the results from the usual plasma model (column 3) deviate significantly 
from the predictions of the generalized plasma model (column 2). At the same time 
the results in column 2 are in close agreement with computations at T = (column 4) 
which are both consistent with experiment. 

Direct comparison with the experimental data of [IE] confirms that the generalized 
plasma model (pQ), (jSJ) is consistent with measurements of the Casimir force at short 
separations. In Fig. 2(a) we plot the differences between the theoretical forces computed 
using the generalized plasma model, F* h , and the mean values of the measured forces. As 
is seen in Fig. 2(a), almost all points are inside the error bars for the difference between 
theory and experiment [12, 28J computed at a 95% confidence level (solid lines). Notice 
that the comparison of data with the zero-temperature theoretical force practically 
coincides with that shown in Fig. 2(a). By contrast in Fig. 2(6) the theoretical force, F^ h , 
is computed using the usual plasma model. It is clearly seen that at short separations 
the usual plasma model is inconsistent with data whereas at larger separations both 
models are in agreement with measurements. 

Refering back to the Introduction, the usual plasma model is consistent with the 
data of experiment p2] which is the most precise experiment at separations from 160 to 
750 nm. At such large separations the predictions of the generalized and usual plasma 
models almost coincide (in fact, the use of the generalized plasma model instead of the 
usual one leads to slightly better agreement between experiment and theory in |12j). 
Thus we can conclude that the generalized plasma model (T5]) both exactly satisfies 
the Kramers-Kronig relations and is also the only model consistent with all precise 
experiments on the Casimir force performed to date. 

4. Discussion 

As was mentioned in the Introduction, the usual plasma model neglects the role of 
free electron scattering. In the absence of scattering the conductivity of free electrons is 
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purely imaginary and the dielectric permittivity is completely real signifying the absence 
of dissipation. The generalized plasma-like dielectric permittivity (JTJ, ([2]), as well as the 
usual plasma model, does not take into account the scattering processes of the free 
electrons. However, it includes dissipation due to the interband transitions of core 
electrons. These transitions are described by a set of oscillators with nonzero resonant 
frequencies. Our results demonstrate that the generalized plasma model exactly satisfies 
the Kramers-Kronig relations and is also consistent with all available experimental data. 
In the same way, as in (111 EU |30], it can also be shown that this model is in agreement 
with the Nernst heat theorem. 

The dissipation processes of free electrons on phonons, impurities etc. are not 
included in the generalized plasma model. They can be taken into account by modifying 
Eq. (DO) to the form 

e(u>) - 1 = A(u) - /°. (18) 
u{u + ig ) 

where /o = oc. The second term on the right-hand side of Eq. (JTBl is the contribution 
of an oscillator with zero resonant frequency, ujq = 0, which was not included in Eq. (j2J). 
This additional oscillator results in a first order pole and leads to the Drude-type term in 
the dielectric permittivity. The Kramers-Kronig relations for the dielectric permittivity 
(fl8~|) are familiar [2]. However, as was discussed in the Introduction, the use of the 
dielectric permittivity (|T8l) leads to a violation of the Nernst heat theorem for perfect 
crystal lattices and to a contradiction with experiment [12]. The question why the 
inclusion of one type of dissipation (interband transitions of core electrons) in the Lifshitz 
theory is necessary while that of another (scattering processes of free electrons) leads 
to contradictions with fundamental physical principles and experiment remains open. 
Future theoretical and experimental developments will shed light on this issue. 
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